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Abstract 

Global characteristics in the diurnal components of OGO-6 neutral mass 
spectrometer measurements near 450 km are discussed qualitatively as well as 
quantitatively on the basis of a theoretical model. Observations and conclusions 
are summarized: (a) During equinox the temperature maximum occurs after 
1600 LT at the equator and shifts toward 1500 LT at the poles, while the oxygen 
concentration at 450 Ian peaks about one hour earlier. This observed phase 
difference is attributed primarily to wind induced diffusion which is most effective 
for the mass density near 250 km and thus can contribute considerably to the 
temperature-density phase anomaly in the thermosphere, (b) there is general 
agreement between the magnitudes and phases of the diurnal, semidiurnal and 
terdiurnal temperature components at 450 km from theox’y as well as OGO-6 and 
radar backscatter measurements, (c) the maximum in the diurnal variation of He 
is observed near 1030 LT consistent with theoretical results which further em- 
phasize the importance of dynamics and diffusion, (d) during solstice conditions 
the diurnal temperature maximum shifts toward later local times, in substantial 
agreement with radar temperature measurements, (e) the temperature-oxygen 
density phase difference at 450 km is observed to decrease with latitude from the 
winter toward the summer hemisphere, where oxygen may even peak after the 
temperature at high latitudes. This effect is attributed to the diurnal variation 
in the O 2 dissociation which becomes most significant in the summer hemisphere 
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where the dissociation rate is enhanced and where the atomic oxygen abundance 
in the lower thermosphere is significantly reduced due to the annual component 
in the large scale circulation. 
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I. INTRODUCTION 


The diurnal variations in the thermospheric density and temperature became 
a new focus of interest when radar backscatter observations revealed a temper- 
ature maximum around 1600 LT (Carru et al. (1967X Mahajan (1969)) which was 
substantially later than the density maximum near 1400 LT as inferred from 
satellite drag data. This discrepancy, commonly referred to as the "phase 
anomaly," constitutes a geophysical problem, whose explanation has been attempted 
by a number of authors ((Rishbeth (1969), Chandra and Stubbe (1970), Mayr and 
Volland (1972a, 1973) and Mayr et al. (1973)). A consensus on the origin of this 
effect has not been reached, primarily because accurate measurements of the 
diurnal variations of temperature and density in the lower thermosphere are not 
yet available to permit convergence on particular mechanisms. 

However, the recent OGO-6 neutral composition measurements (Carignan 
and Pinkus, 1968) near 450 km complement radar backscatter and satellite drag 
data and provide important additional inputs to the phase anomaly problem. 
Presumably, N 2 is closest to diffusive equilibrium. Its variation is thus con- 
sidered to be a good measure of the variation in the gas temperature although 
this assumption will be qualified further in this paper. Differences in the time 
response of various constituents and their relationships to the temperature (Nj) 
phase can hence be inferred. Furthermore, the good data coverage in local time, 
latitude and season permits an empirical description of these data with relatively 
high temporal and spatial resolution (Hedin et al. 1973a, b). 
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In this paper we shall discuss the local time variations inferred from the 
OGO -6 neutral mass spectrometer data. The seasonal components will be 
treated separately (Reber et al. 1973). 

II. OGO -6 NEUTRAL COMPOSITION MEASUREMENTS 

1. Empirical Model 

Carignan and Pinkus (1968) described the neutral mass spectrometer and 
its operation on OGO- 6 . Global properties in the distribution of composition 
and temperature, inferred from these measurements, were reported by 
Reber et al. (1971), Taeusch et al. (1971), Hedin and Reber (1972), Carignan 
and Reber (1971), Hedin et al. (1972, 1973a, 1973b), Hedin and Mayr (1973). 

A detailed account of the statistical analysis for the empirical model of 
the OGO -6 mass spectrometer data is given in a separate paper (Hedin et al., 
1973b). Our discussion will be based on the local time components in this 
OGO -6 model; thus a brief review of its properties is given here: 

(1) The composition data for the model are sampled for magnetically quiet 
conditions (Ap 4 7, ap 12) in the 400-600 km altitude range. 

(2) Universal time and longitudinal effects are eliminated by considering 
only longitudinal averages. 

(3) The observed N variations are interpreted in terms of T (temperature) 
variations whereby the assumptions are made that N, and T„ are invariant 

at 120 km and that the shape of the height distribution of T^ is known, with T^ 
(exospheric temperature) being a free parameter to be determined from the 
observed N 2 concentration. 
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(4) Assuming diffusive equilibrium for all constituents at altitudes above 
120 km, the inferred temperature distribution is adopted to describe the composi- 
tion of N^, O and He at 450 km (representative height of the measurements) and 
to deduce variations in O and He at 120 km (assuming N 2 is invariant there). 

(5) The atmospheric parameters , O and He at 450 km and O and He 

at 120 km are represented in an expansion of spherical functions including diurnal, 
semidiurnal and teridiurnal harmonics for the local time variations, and annual 
and semiannual harmonics for the seasonal variations. 

The model composition at 450 km can be regarded as a good representation 
of the OGO-6 measurements. Density extrapolations to this altitude level cover 
only about one scale height for each species and thus any uncertainties in the 
temperature determination are not very significant. However, two characteristics 
of the empirical model must be kept in mind: a. the temperatures are not 
measured directly but are inferred under assumptions (see point three of the 
previous paragraph) which are not completely justified on theoretical grounds 
(Volland and Mayr (1970, 1973), Mayr and Volland (1973), Mayr et al. (1973)) and 
b. the assumption of diffusive equilibrium is not completely valid below 200 km 
(see also Mayr and Volland (1971, 1972b), Reber and Hays (1973)) thus the 
inferred changes in the neutral composition at 120 km are probably unrealistic 
particularly with regard to the diurnal components which cannot be easily excited 
over such relatively short time periods. For these reasons, in particular, the 
OGO-model cannot be a substitute for in situ temperature and composition 

measurements in the lower thermosphere. 
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On the other hand, theoretical investigations have shown that no matter how 
complex the dynamic properties and the phase relations of composition and 
temperature are in the lower thermosphere, at exospheric heights where the 
OGO-6 data are sampled, the temperature and the N 2 concentration vary nearly 
in phase (Mayr et al., 1973). Thus the observed N 2 variations at 450 km are in 
practice a good measure for the variations in the exospheric temperature, 
and consequently the inferred composition should also be reasonable at least 
down to 200 km where diffusive equilibrium holds. In our discussion of the 
diurnal variations we shall therefore restrict ourselves primarily to altitudes 
above 200 km which, incidentally, is also the height region in which most satellite 
drag data have been taken. 

2. Observations 

Figures 1 through 4 show contour plots for the local time components of , 
N 2 , O and He at 450 km for September equinox and June solstice conditions taken 
from the OGO-6 model (Hedin et al., 1973b). The temperature plots describe the 
deviations from the global mean (1098 °K). The density plots represent the natural 
logarithms of the ratios of the variations to the global means (2.16 x 10® for N 2 , 
6.72 X 10^ for O and 3.02 x 10® for He). In reality, seasonal and time independent 
latitudinal variations have to be considered superimposed on these local time com- 
ponents which of course, would not change the phase distributions for the 
various atmospheric parameters. The latitudinal variations in the times of 
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density maxima are indicated with solid lines. For comparison the distributions 
for the temperature phase are shown in dashed lines. The diurnal variations in 
the OGO-6 model are not well defined during solstice conditions at high latitudes 
and areas of poor data coverage are shaded in Figures 1 through 4. 

Considering the uncertainties in the coefficients of the empirical model we 
estimate the errors for the magnitudes of the diurnal variations as ±15% and for 
the times of maxima as ±0.5 hr. This phase imcertainty is small compared with 
the observed phase differences between He and T^ . However, the phase difference 
between O and T^ is only marginally determined. 

HI. QUALITATIVE DISCUSSION 

Figures 1 through 4 reveal a number of interesting features in the global 
distribution of the thermospheric composition. In the following, some of the 
outstanding effects are classified and discussed: 

1. "Equatorial anomaly" 

During conditions of equinox it is apparent that the day time bulge and the 
night time minimum in the density and temperature are elongated around 
the equator. Shallow day time maxima and night time minima are observed 
near ±20° latitude indicating maximum diurnal amplitudes at these latitudes. 

This phenomenon resembles the equatorial anomaly in the Fj region (see the 
review paper by Burkard,1966) which has been attributed to ExB drifts (Hanson 
and Moffett, 1966) induced by electric fields from the dynamo region of the 
lower thermosphere. 
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The connection with the equatorial anomaly has been demonstrated in 


Hedin and Mayr (1973). It was shown there that (a) the shallow equatorial day 
time minimum and night time maximum in the gas temperature tend to follow 
the magnetic equator when individual satellite orbits are considered over various 
longitudes and (b) the variations of ion drag associated with the equatorial anomaly 
can affect the thermospheric circulation and energy advection sufficiently to 
produce the observed structure in the gas temperature. From Figure 3 it is 
apparent that the atomic oxygen bulge is also elongated around the equator thus 
following to some extent the temperature pattern. 

He in Figure 4 on the other hand shows a pronounced day time peak at the 
equator and two minima at 0200 and around 2100 LT during night. This is in part 
a consequence of the semidiurnal and terdiurnal components which are most 
important at low latitudes although the errors in determining these modes are 
relatively large (30%). Nevertheless, this tendency appears to be consistent 
with our picture in which the distribution of this species is significantly controlled 
by diffusion effects induced by the thermospheric circulation (Mayr and 
Volland, 1973). This process is affected by ion drag in that any decrease in 
the electron density at the equator would increase the wind velocity. This in 
turn would enhance the diffusion process, which, if dominant, would explain 
then the increase in the diurnal amplitude of He. 
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2. Phase Differences and Their Latitude Dependence. 

Figures 1 through 4 reveal a strong latitude dependence and systematic 

differences between the phases of the indiAddual constituents. At low latitudes 
and during equinox conditions it is apparent that (and ) peak latest in the 
afternoon (1600 LT), preceded by O and He, the latter reaching its maximum 
already before noon. The magnitude of this effect is largest in the equatorial 
region where the T -O and the T -He phase differences are about one and six 

e g 

hours respectively. At high latitudes, however, the temperature varies more 
or less in phase with atomic oxygen, reaching its peak earlier in the afternoon 
at 1400 LT, while the He maximum shifts toward later local times (1200 LT) at 
high latitudes thus reducing there T^ -He phase differences to two hours. Although 
there is a half hour uncertainty in the determination of the diurnal density maxima, 
conceivably affecting some details in the global distribution, the observed trends 
are probably real. 

Two factors seem to contribute. One is related to the global distribution 
of the solar heat input which affects directly the temperature phase, and the 
second one is related to the global thermospheric circulation which influences, 
through diffusion, primarily the lighter constituents. 

The diurnal variation in the energy input, which can be treated as an 
impulse with a half width of 2-5 hours, produces a similar but skewed temperature 
variation with the peak occurring close to the time where the heat input starts 
to decrease. At the equator, where absorption effects are least significant, the 
heat input is relatively broad around noon and hence the temperature peaks 
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late in the afternoon. Toward higher latitudes absorption effects become in- 
creasingly import (thus narrowing the effective width of the heat source) and there- 
by shift the temperature maximum toward earlier local times in agreement with 
the phase distributions in Figures 1 and 2. In the OGO-6 model this property 
is essentially described with the higher harmonics which are most significant 
at low latitudes. There in particular, the terdiurnal component plays the role 
of shifting the temperature peak toward the late afternoon. The reason for this 
is that the solar heat input, with its flat and elongated maximum around noon, 
has a relatively strong third harmonic which peaks at 1600 LT, in contrast to 
the energy input maxima near noon in the first two harmonics. 

Wind induced diffusion causes redistribution of the minor constituents in the 
lower thermosphere. In this process the lighter constituents O and He are 
transported from late afternoon toward early morning hours such that their diurnal 
peaks are shifted away from the temperature maximum toward earlier local times 
(Mayr and Volland, 1972, 1973a). This effect is more important for He than 
for O. Furthermore, the temperature induced variation in the composition, which 
competes with the diffusion effect and which tends to produce density variations 
in phase with the temperature, is more effective for the heavier constituent O. 

For both these reasons the He-Tg (N 2 ) phase difference should be larger than 
that between O and Tg in substantial agreement with the OGO-6 observations 
in Figures 1 through 4. 
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The observed latitude dependence in the phase differences between O, He 
and Tg are consistent with certain dynamic properties of energy advection and 
diffusion which are related to the horizontal scales in the global thermospheric 
circulation. Diffusive mass transport in O and He is to a great extent proportional 
to the vertical wind velocity. This in turn increases with increasing divergence of 
the horizontal wind field. For a given horizontal wind velocity the diffusion 
effect thus increases inversely proportional to the scale of the circulation or 
with the number of circulation cells (see Mayr and Volland, 1973b). At low 
latitudes, where the higher harmonics in the diurnal variations are most 
significant, the diffusion process is thus most effective, causing the phase 
differences between He, O, and Tg to increase toward the equator (disregarding 
the low latitude structures induced by the equatorial anomaly in the F 2 - region). 

In parallel, heat advection, with global characteristics similar to diffusion, 
also becomes more effective toward the equator thus damping there the tempera- 
ture amplitude relative to what it would be in the absence of this energy transport. 
It follows then that the temperature induced variations in the composition are 
relatively less significant at low latitudes hence further contributing there to the 
increased phase difference between the gas temperature and the lighter 
constituents (Mayr et al. 1973). 

3. Seasonal Effects. 

Figures 1 through 3 show that the diurnal variations for O and Tg (N 2 ) are 
significantly different in the summer and winter hemispheres, while the He 
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distribution in Figure 4 is fairly symmetrical with respect to the equator. We 
recognize that the atmosphere is poorly defined by the OGO-6 model above 50° 
latitudes. Nevertheless, even at low to mid-latitudes there are a number of pro- 
nounced features in the composition such as; (1) the increasing Tg (1^) amplitude 
between +30° and -30° latitudes from the summer toward the winter hemisphere; 

(2) the broadenir^ of the diurnal Tg (N^) and O maxima in the summer hemisphere 
accompanied by a shift of the Tg (Nj) peaks toward later local times near the 
subsolar point; and (3) the decrease of the phase difference between O and Tg in 
the summer hemisphere with some indication that oxygen might even peak later 
than the gas temperature at high latitudes . 

A number of factors contribute to these effects, some of which we shall 
further discuss quantitatively in the subsequent chapter of this paper. In the 
Slimmer hemisphere the time span of solar illumination increases thus producing 
there the observed broadening in the diurnal maxima of oxygen and temperature 
while shifting their peaks toward later local times. With the widening of the 
solar heat input, the ratio of the diurnal to the local time average components 
decreases in summer relative to winter thus contributing to the winter maximum 
at -30° latitude in the diurnal temperature amplitude. The actual magni- 
tudes of these effects are probably also affected by the seasonal modulation 
of the diurnal wind circulation as well as by the seasonal (winter to summer) 
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variations in the local time average components of temperature and composition. 
While the former tends to damp the horizontal variations in Tg, the latter may 
substantially contribute to the observed increase of the temperature amplitude 
in winter when the average background pressure (internal energy) is lower 
relative to the stunmer hemisphere. 

The seasonal component in the diurnal variations of atomic oxygen, and in 
particular its phase relation to Tg, is governed by a complex superposition 
of transport processes related to temperature variations, wind circulation and 
dissociation of Oj 

O 2 + hv - 0 + 0 

Among these processes, dissociation of probably contributes most to the 
seasonal variations. The significance of this reaction is proportional to the 
dissociation rate 

Qms'"2F[Oj (2) 

(where Fis the dissociating solar flux in the Schumann Runge continuum and [O 2 ] is the 
concentration of O^) and it is inversely proportional to the amount (reservoir) of 
atomic oxygen. From winter to summer the dissociation rate increases while the 
local time average atomic ojsygen concentration decreases owing to mass transport 
associated with the annual and semiannual components in the global circulation (Mayr 
andVolland, 1971, 1972b). For both these reasons the dissociation process should 
thus be much more important for the diurnal variations in summer than in 
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winter. With dissociation maximizing around noon its cumulative effect tends 
to produce a diurnal maximum in O during the late afternoon, thereby shifting 
the oxygen phase away from the temperature toward later local times. This 
trend is exactly opposite to the wind induced diffusion process, discussed earlier, 
which tends to shift the oxygen phase away from the temperature toward earlier 
local times, and both processes are different from the temperature induced 
variations which would tend to keep atomic oxygen in phase with the gas 
temperature. 

IV. QUANTITATIVE DISCUSSION 

Our quantitative discussion will be based on a quasi -three dimensional 
model (Mayr et al. 1973). Applying perturbation theory the equations of energy, 
mass and momentum conservation are solved for a two component gas whereby 
the momentum transfer between the two constituents (diffusion) is considered 
in a self consistent form. In particular two versions are being used: one in 
which oxygen diffuses through Nj (+O 2 ); and another one in which He diffuses 
through a fictitious gas component with mass and density corresponding to the 
mixture of O, N^ and . For He, thermal diffusion has also been considered. 

We consider the first version to be realistic with regard to temperature and 
wind fields since the energetics of the thermosphere is essentially determined 
below 500 km where O and N^ are the major constituents. The second version, 
however, serves only to give an estimate for the diurnal variations of He and 
is probably not realistic energetically. 
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From a Fourier analysis of solar heat input calculations by Harris and 
Priester (1962) the relative magnitudes and height distributions of the diurnal, 
semidiurnal and terdiurnal components are adopted as shown for the equator in 
Figure 5. The absolute magnitudes are determined such that the resulting theo- 
retical temperature amplitude in the first diurnal component agrees with OGO-6 

observations. This model of the heat input distribution deviates from the one 
in Mayr et al. (1973) in that the amplitude of the third harmonic is considered 
constant below 200 km reflecting the importance of absorption. We retain, however, 
all other assumptions regarding the latitude dependence and the phase variations 
with height. Furthermore, we have extended the earlier analysis by including 
dissociation (Equation 1), Figure 5 shows the amplitudes of the dissociation 
rates for the diurnal and semidiurnal components at the equator. Considering 
that dissociation is only significant in the lower thermosphere and mesosphere 
we neglected its third harmonic consistent with the low heat input rates in this 
component. It is assumed that the maxima in the diurnal components of the 
dissociation rate occur at 1200 LT and that the amplitudes decrease with 
latitude proportional to sin 6 and sin^^ for the diurnal and semidiurnal com- 
ponents respectively (0 = colatitude) . 

For the local time average ion drag an electron density distribution was 
adopted with an F 2 maximum of 1.3 x 10®/cc at 300 km and an electron density 
of 3.5 X 10"^ at 120 km. The diurnal component of the electron density was 
assumed to vary with a relative amplitude of 0.5 and a maximum at 1200 LT. 
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All other input parameters such as the coefficients for heat conduction, viscosity 
and molecular and eddy diffusion were the same as in previous papers (Mayr 
and VoUand, 1973aj Mayr et al. 1973). 

1. Tidal Components 

The computed diurnal variations in Tg , mass density and He concentration 
are presented in Figure 6 for the equator. For comparison are shown the 
OGO-6 Tg, oxygen (mass density) and He measurements at 450 km, as well as 
the diurnal variations in the mass density inferred from satellite drag 
(Jacchia, 1971) under similar solar activity conditions. To simplify the presen- 
tation the theoretical He distribution is only given for 450 km while the mass 
density and temperature variations are shown for various heights. 

The agreement in the magnitude of the diurnal temperature variation is 
not significant since the heat input for the first diurnal component was adjusted 
such that the amplitude (not phase) in the theoretical temperature component 
matched the corresponding Tg tide in the OGO-6 measurements. However, 
the relative amplitudes of the diurnal, semidiurnal and terdiurnal tides as well 
as the temperature phase and its relation to the phases of He and O at 450 km 
are significant, since in the framework of the theoretical model (Mayr et al. 1973) 
no adjustments or assumptions have been made to control these properties. 

From Figure 6 it is apparent that the maximum in the theoretical temperature 
variation occurs at 1600 LT in close agreement with the OGO-6 observations. 
Furthermore the computed maximum in the mass density near 1500 LT is in 
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substantial agreement with both the OGO-6 and satellite drag data. For He the 
phases agree reasonably well in particular for the dominant first harmonic which 
produces a maximum near 1100 LT in the theory. Due to the second and third 
harmonics, however, an elongated He bulge develops which has a shallow maximum 
in the leading edge after 1200 LT. However, in view of relatively large errors 
associated with the semidiurnal and terdiumal components of the empirical He 
model from OGrO-6 the discrepancy between theory and observation may not 
be significant. 

The physical processes for the phase shift of He and oxygen away from the 
temperature toward earlier local times, revealed in the theoretical results in 
Figure 6, have been quantitatively discussed in some detail (Mayr and Volland 
(1972, 1973a), Mayr et al.(1973)). In this theory it is energy and mass transport 
associated with the thermospheric circulation which contribute to the phase 
differences between individual constituents by (a) producing a phase shift in the 
temperature toward earlier local times at lower altitudes and (b) inducing devia- 
tions from diffusive equilibrium in the lower thermosphere such that the lighter 
and minor constituents tend to peak 90° (6 hours for the diurnal component) 
before the temperature maximum. 

This is illustrated in Figure 7 where the theoretical phase distributions 

are shown together with the phases of mass density and temperature deduced 

from the OGrO-6 measurements. The temperature determination from Nj 

measurements on OGO-6 is based on the assumption that Nj is invariant at 120 km 

and that the temperature phase is height independent. Thus the N 2 observations 

at 450 km can only produce an "effective temperature phase" representative 
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of the entire thermospheric temperature structure, even if Nj is in perfect 
diffusive equilibrium. In reality, the Nj maximum may occur either before or 
after the Tg peak at exospheric heights depending on whether the temperature 
phase shifts at lower heights toward earlier or later local times respectively. 

On the basis of the theoretical model discussed here the temperature should 
therefore actually peak after the density. The N^-Tg phase difference, how- 
ever, is small due to (1) the large influence of the temperature phase above 200 Ian 
which is almost constant and (2) due to the drag (diffusive) interaction with O which 
tends to shift the phase toward later local times, thereby partially compen- 
sating the effect of the earlier temperature phase at lower heights. With this 
qualification in mind, the extrapolation of OGO-6 temperature and density varia- 
tions to lower heights is not very meaningful below 200 km, particularly not for 
O and He which are not in complete diffusive equilibrimn there. Down to this 
altitude theory and "observations" agree within the ±0.5 hour uncertainty in the 
phase determination at 450 km. The agreement is particularly good in the 
magnitude of the temperature-density phase difference as well as its height 
dependence, which clearly indicates a widening of the phase discrepancy at lower 
altitudes. Since it is inherent in the present analysis of satellite drag data that 
the density phase is height independent it is then likely that these observations 
only produce a phase representative of some intermediate altitude region around 
300 km, where we would also put the density phase at 1430 LT. 
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The density extrapolation from the OGO-6 measurements to greater heights 
suggests that the density-temperature phase difference decreases further up to 
600 km where it starts again to widen due to the growing importance of He. This 
reversion in the phase distribution is particularly significant in the winter 
hemisphere where the He content is considerably enhanced; likewise it is least 
important in summer where the He density reaches its minimum. 

Magnitudes and phases of the three diurnal temperature harmonics are 
compared in Table 1. Considering that in our theory and P^ are the 

spherical harmonics associated with the diurnal, semidiurnal and terdiurnal 
components, we have chosen only the corresponding coefficients from the OGO-6 
model. Similarly, the radar temperature measurements reported by Salah and 
Evans (1973) at 45® were extrapolated to the equator adopting the latitude de- 
pendences P J (siad), P2 (sin2 6) and P| (sirf<9) for the three harmonics. 

Table 1 reveals generally good agreement between the three data sets. The 
exception is the magnitude of the semidiurnal tide, inferred from OGO-6, which 
is by about a factor of two lower than those inferred from theory or radar 
temperature measurements. Although the errors in the OGO-6 measurements 

of the semidiurnal tide are relatively large (20%) they cannot completely account 
for this disagreement. This discrepancy may be due to the fact that the OGO-6 
temperature model is not based on temperature but on Nj measurements. On 
the other hand the theoretical description of the semidiurnal component is 
problematic in that nonlinear effects and energy coupling from the lower atmosphere 
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play a major role. Considering the extrapolation procedure adopted for the 
radar temperature components their agreement with the theory is perhaps not 
convincing enough either to rule against the small semidiurnal component in the 
OGO-6 data. 

2. Seasonal Variations 

Figure 8A shows the latitudinal variations in the local times of maxima 
for , O and He from the OGO-6 observations during June solstice and September 
equinox. Figure 8B shows the phase differences for O and He with respect to T^. 
While for He there is relatively little change between summer and winter, the 
-O phase difference reverses from the winter toward the summer hemisphere. 
This feature is further characterized in Figure 8C where the difference in the 
times of maxima between equinox and solstice conditions is shown for O. From 
this it is apparent that relative to equinox the ojQ^gen phase shifts toward later local 
times in summer and toward earlier local times in winter. 

We suggest that this feature is caused by the diurnal variation in O 2 dissocia- 
tion as it is seasonally modulated by the winter to summer decrease of the oxygen 
content of the lower thermosphere (see the qualitative discussion in this paper). 

To demonstrate the significance of this effect a parametric study has been 
carried out for mid-latitudes, whereby the diurnal amplitude of the dissociation 
rate, as well as the oxygen concentration at 120 km have been varied independently. 
Figure 9 shows the diurnal component of theTg-O phase difference at 450 km 
as a fimction of these two parameters. 
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For zero dissociation rate the phase difference is almost an hour at 450 km 
(and almost two hours at 250 km). As has been pointed out before, this phase 
difference is to a great extent generated by diffusion and for this process it is 
important that oxygen is lighter than the mean molecular mass and is a minor 
constituent in the lower thermosphere. With increasing oxygen concentration the 
latter condition deteriorates, thus diminishing the importance of diffusion. Hence 
the Tg -O phase difference slightly decreases, as can be seen from Figure 9 
(when dissociation is neglected). 

The diurnal variations in Oj dissociation, which produce atomic oxygen 
at a rate that peaks during noon, induces diurnal variations inO which tend to 
peak near dusk. With the production rate being independent of 0,the significance 
of dissociation increases as the background density (local time average) of 
atomic oxygen decreases. The oxygen phase is thus delayed toward later local 
times when, for a particular dissociation rate, the oxygen concentration is decreased 
and or when the dissociation rate is increased for a given oxygen concentration. 

Following an interpretation of the F 2 -region winter anomaly by Johnson (1964) 
and King (1964) it has been shown theoretically that the annual component in the 
thermospheric wind field induces about a factor of three decrease in the atomic 
oxygen concentration of the lower thermosphere between winter and summer 
(Mayr and Volland (1971, 1972b)). This seasonal effect in the neutral composition 
is found to be consistent with the OGO-6 neutral mass spectrometer measurements 
at 450 km (Hedin et al. 1973, Reber et al., 1973). Figure 9 should therefore be 
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read in the sense that lower oxygen concentrations (and higher dissociation rates) 
are representative for summer conditions while higher oxygen concentrations (and 
lower dissociation rates) represent winter conditions. For a wide range of dissocia- 
tion rates (between 2 x 10“®and 5 x lO'^at 120 km) a factor of three decrease in 
atomic oxygen thus changes the Tg-O phase difference by about an hour, from oxygen 
peaking substantially before the temperature in winter toward oxygen being in 
phase or peaking slightly after the temperature in summer. Trend and magnitude 
of this effect are thus in substantial agreement with the observed phase distribu- 
tion for O revealed in Figure 8. 

These results are summarized in Figure 10 for which we computed the diurnal 
variations of temperature and oxygen concentration at 450 km in the summer and 
winter hemispheres. Inputs on the local time averages of temperature and 
composition were taken from the OGO-6 measurements (Hedin et al. 1973b) which 
show that at 120 km O should (a) decrease by a factor of 1.5 from the equator 
toward the summer hemisphere and (b) increase by a factor of 1.3 from the 
equator toward the winter hemispheres at midlatitudes. The exospheric tempera- 
ture is observed to change from 1098°K at the equator to 1170°K and 910°K in 
the summer and winter hemisphere respectively. The density variations in 
Figure 10 are normalized to 1 for the local time average equatorial oxygen con- 
centration at 450 km. For the heat input rate it was assumed that the three 
harmonics vary with latitude like sin {6- 6J, sin^{9 - 6^ ), sin^ (9 - 9^), 
where 9^ is the colatitude of the solar declination. The latter assumption was 
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made to account for the fact that absorption effects are least significant 


at subsolar latitudes and thus higher harmonics, particularly the terdiurnal 
component, are most significant there. It is apparent from Figure 10 that the 
theory can basically reproduce the phase shift of Tg toward later local times in 
summer and the narrowing of the density-temperature phase difference from 
winter to summer; however it cannot reproduce the large diurnal amplitude in 
winter, a feature which may in part be related to the winter enhancement in the 
F 2 region ionization (winter anomaly). 


V. SUMMARY AND CONCLUSION 

Global characteristics in the diurnal components of OGO-6 neutral mass 
spectrometer measurements (Carignan and Pinkus, 1968) near 450 km, which 
have been described empirically in analytic form (Hedin et al. 1973b), are 
discussed qualitatively as well as quantitatively on the basis of a theoretical 
model (Mayr et al. 1973). 

The observations and conclusions are summarized: 

(a) During equinox the diurnal temperature maximum occurs after 1600 LT 
at the equator and shifts toward 1500 LT near the poles due to the diminishing 
importance of higher harmonics (particularly the terdiurnal one) which describe 
the broadening of the diurnal solar heat input at subsolar latitudes and due to 
the decreasing ion drag. 
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(b) During equinox the oxygen concentration at 450 km is observed to peak 
near 1500 LT at the equator and shift toward 1430 LT at the poles. This phase 
difference between and O is consistent with a theoretical prediction (Mayr 
and Volland, 1972) that attributed this effect primarily to wind induced diffusion. 
This Tg -O phase difference increases toward lower altitudes, thus producing a 
phase difference between T^ and mass density of about 1.5 hours at 250 km 
which contributes considerably to the phase anomaly in thermosphere. 

(c) From observations as well as theoretical analysis it is concluded 
that the phase of the mass density as well as its difference with respect to the 
temperature can vary by more than an hour between 200 and 600 km and by 
almost an hour between equator and poles during equinox conditions. 

(d) There is general agreement between the magnitudes and phases of 
the diurnal, semidiurnal and terdiurnal temperature components at 450 km 
from theory as well as OGO-6 and radar backscatter measurements, except 
for the magnitude of the semidiurnal tide which is somewhat small (2.1% 
at the equator) in the OGO-6 analysis. In particular, the observed large 
amplitude (3%) and late phase (1630 LT) in the terdiurnal tide are reproduced 
theoretically suggesting that it is excited directly by the solar UV component 

associated with the elongation of the heat input maximum around noon. 

(e) The maximum in the diurnal variation of He is observed near 1030 LT, 

in substantial agreement with the theoretical results, thus emphasizing the 
importance of wind induced diffusion which primarily accounts for this large 
phase difference with respect to Tg. 
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(f) From winter to summer the temperature maximum is obseiwed to shift 
toward later local times in ^reement with radar temperature measurements 
(Salah and Evans, 1973). This effect is attributed to the broadening of the solar 
heat input maximum in summer, particularly near subsolar latitudes where 

the time period of illumination is long and absorption effects are least significant. 

(g) Dviring solstice conditions it is observed that the positive phase difference 
between Tg and O decreases from the winter toward the summer hemisphere 
with some indication that O might even peak after Tg at high latitudes. This 
effect, which indicates significant variations in the phase of the mass density 
between winter and summer, is attributed to the diurnal variations in O 2 
dissociation which produces oxygen during daytime and thus tends to shift 

the O peak toward later local times. Between winter and summer the significance 
of dissociation increases due to the increasing dissociation rate but most 
importantly due to the decreasing O concentration in the lower thermosphere. 

The greater importance of dissociation in summer thus compensates (or even 

overcompensates) the wind induced diffusion which normally shifts the oxygen 
peak away from the temperature toward earlier local times. 

Important insight into the thermosphere structure has been gained from the 
OGO-6 composition measurements around 450 km. However, to converge with 
some degree of uniqueness upon the mechanisms responsible for these 
structures it will be necessary that similar global observations, combined 
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with in situ temperature and wind measurements, be carried out at lower 
altitudes where the dynamics and energetics of the thermosphere, as well as 
coupling with the mesosphere, have the greatest impact. 
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LOCAL TIME 


Figure I. Diurnal variations in the exospheric temperature inferred from the OGO-6 N2 measurements around 450 km. The 
contour plots represent temperature differences from the local time average component for September equinox (left) and June 
solstice (right) conditions. 
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Figure 2. Contour plots of the logarithm (base e) of the ratio of the 450 km N2 densities to the local time average component 

for equinox and solstice conditions. 




Figure 3. Contour plots of the logarithm (base e) of 
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ratio of the 450 km O densities to the local time average component 
and solstice conditions. 
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Figure 4. Contour plots of the logarithm (base e ) of the ratio of the 450 km He densities to the local time average component 

for equinox and solstice conditions. 
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Figure 5. 


Height distribution of the amplitudes of the diurnal components of the heat input Q 
(solid lines) and dissociation rates (dashed lines) D at the equator. 
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Figure 6. The theoretical diurnal variations in He at 450 km and in mass density and temperature 
at various heights are shown for the equator during equinox. At 450 km these results are com- 
pared with the corresponding OGO-6 oxygen and satellite drag mass density measurements. 
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Figure 7. Height distribution of the phases of mass density and temperature from theory and 
inferred from OGO-6 neutral composition measurements around 450 km. Below 250 km the 
OGO-6 extrapolations become uncertain as indicated with dashed lines. 
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Figure 8. (a) Latitudinal distribution in the local times of maxima for 0 and He (heavy lines) and 
Tg (thin lines) at June solstice (solid lines) and September equinox (dashed lines) from the OGO-6 
measurements at 450 km. (b) T^-O and "^^a’lde phase differences for June solstice (solid) and 
equinox (dashed) conditions, (c) Phase difference between June solstice and September equinox 
conditions of 0. 
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Figure 9. Phase difference for the diurnal component at 450 km as a function of the diurnal 
amplitude in the dissociation rate (ordinate) and the atomic oxygen concentration at 120 km. 
Note the change in the phase sequence for low oxygen concentrations and high dissociation 
rates corresponding to summer conditions. 
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Figure 10. Diurnal variations in and 0 for summer (solid lines) and winter (dashed lines) hemispheres from theory and 
OGO-6 observations at 45° latitudes. The local time average temperatures and oxygen concentrations were adopted 
from the OGO-6 model (Hedin et al., 1973). 
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